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We present an investigation of the inclusive semileptonic decay widths of the heavy baryons Aq and Sq (Q — b, c) 
performed within a relativistic constituent quark model, formulated on the light-front. In a way conceptually similar to the 
deep-inelastic scattering case, the f/Q-baryon inclusive width is expressed as the integral of the free Q-quark partial width 
multiplied by a bound-state factor related to the Q-quark distribution function in the Hq. The non-perturbative meson 
structure is described through the quark-model wave functions, constructed via the Hamiltonian light-front formalism 
using as input the Kalman-Tran-D'Souza equal time wave functions. A link between spectroscopic quark models and the 
Hq decay physics is obtained in this way. It is shown that the bound-state effects and the Fermi motion of the 6-quark 
remarkably reduce the decay rate with respect to the free-quark result. Our predictions for the BR(A C — > X s £ve) and 
BR(A b — ■> X c lvf) decays are in good agreement with existing data. 
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1. Introduction 

The investigation of inclusive semileptonic Aj, decays 
can provide relevant information on the Cabibbo- 
Kobayashi-Maskawa (CKM) parameters \V c b\ and 
\V u b I as well as on the internal non-perturbative struc- 
ture of the Ah baryon. The first evidence of semilep- 
tonic Afc decays had been reported by ALEPH and 
OPAL collaborations who had seen an excess of corre- 
lated A s £~ pairs over A+ pairs from Z decays [fij , 0. 
More recently, the CDF collaboration [|| measured 
the Af, lifetime using its semileptonic decay mode. It 
is therefore important to study theoretical models for 
the semileptonic Aj, decays. 

As far as the theoretical point of view is con- 
cerned, the QCD -based operator product expansion 
(OPE) combined with the heavy quark expansion 
is widely recognized as a consistent dynamical ap- 
proach for investigating inclusive heavy-flavour de- 
cays [Q. Existing theoretical predictions for inclu- 
sive weak decays are in reasonable agreement, within 
the expected range of uncertainty, with the data on 
decays of charmed and beauty particles. However, 
the result of the OPE approach for the A& lifetime 
is puzzling because it predicts that (ta 6 /tb)opb = 
0.98 + 0(1/ 'ml), whereas the experimental finding 
suggest a very much reduced fraction (rA b /TB)exp — 
0.8 ±0.04. It is not clear whether the present contra- 
diction between the theory and the data on t^/tb 
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is a temporary difficulty, or evidence for the funda- 
mental flaws in the OPE approach. In spite of great 
efforts of theoretical activity the A& lifetime remains 
significantly low which continues to spur theoretical 
activity. In this respect, the use of phenomenological 
models, like the constituent quark model, could be 
of great interest as a complementary approach to the 
OPE method. 

What criteria should be used in determin- 
ing an appropriate constituent quark model for the 
baryons? There is now a vast amount of data on 
the spectrum of baryons. Many quark model calcu- 
lations of the spectrum attempt to account for the 
masses, spin etc., without studying the decays. A se- 
rious drawback to this approach is that it is difficult 
to assign predicted baryons to experimentally known 
particles on the basis of spectrum alone. In this sense 
this paper is a continuation of the work of Kalman 
and Tran || to produce a complete picture of the 
spectrum and decays of the baryons. This model is 
based upon the following principles: i) consistency 
with the model of the mass calculation - A change 
of some feature of the mass calculation in the decay 
model may compromise the purpose; ii) The decay 
calculation should be mainly determined by the mass 
calculation. For instance, if there are too many pa- 
rameters in the decay model it is difficult to relate 
the results obtained to the particles considered in the 
model, used in the mass calculation. The decay cal- 
culations were based upon the quark pair creation 
model of Le Yaouanc et. al. || They studied the 
decay processes N, A — > A^7r and A, £ — > NK. The 
pair creation strength 7 is replaced by k 7 . Besides 
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the meson radius, these are the only parameters of 
the model. The results are in good agreement with 
the data. To calculate decay widths in their model, 
Kalman and Tran had to re-correct the originally cal- 
culated parameters a p and a>\ by approximately 0.05 
GeV from the values determined from their spectrum 
calculation, in order to make the decay widths match 
experimental data. The resulting change in a p and 
a\ causes a shift in the originally calculated mass 
spectrum. It is these values that must be used if the 
decay widths are to be correct. 

At the time that the Kalman- Tran model was 
introduced, there was little data available for the b 
and c sectors. The model cannot fit these sectors 
properly without the modifications performed in the 
Kalman and D'Souza model The essential mod- 
ification is that the coupling constant of the interac- 
tion is made to be scale dependent in a manner that 
does not require knowledge of the behaviour at low 
energy scales and does not introduce any new param- 
eter. The calculations of the inclusive semileptonic 
widths and the corresponding branching ratios found 
here introduce no new parameters. It is based solely 
upon the eigenfunctions for baryon states found by 
Kalman and D'Souza |l3| 

In what follows we assume that, instead of 
QCD with its complicated dynamics of infinite num- 
ber of degrees of freedom in the light cloud, we con- 
sider a constituent bound-state problem of a heavy 
quark interacting with a lighter ones via a potential. 
Then, using the formalism of the light-front (LF) rel- 
ativistic quantum mechanics it is possible to encode 
all the nonperturbative QCD effects in a LF quark 
model wave function ip(x, p 2 ) of a heavy hadron. The 
internal motion of a heavy Q-quark inside the heavy 
flavor meson is described by the distribution function 
\ijj(x,p 2 )| 2 , which represents the probability to find a 
heavy quark carrying a LF fraction x=pg /Pjl} Q of the 
meson momentum and a transverse relative momen- 
tum squared p 2 =p 2 . In Ref. this formalism has 
been used to establish a simple quantum mechani- 
cal relation between the inclusive semileptonic decay 
rate of the heavy hadron and that of a free heavy 
quark. In this paper we shall compute the preasymp- 
totic effects for the inclusive semileptonic decays of 
the heavy baryons Hq in the framework of the LF 
quark model, which is a relativistic constituent quark 
model based on the LF formalism. 

2. Basic formulae 

The approach of relies on the idea of du- 
ality in summing over the final hadronic states. It 



has been assumed that the sum over all possible 
charm final states Xq> can be modelled by the de- 
cay width of an on-shell Q quark into an on-shell Q' 
quark folded with the Q-qu&rk distribution function 
f2 Q (x,pf) — \ip < ^ [Q (x,p 2 )\ 2 The latter represents the 
probability to find a Q quark carrying a LF fraction 
x of the hadron momentum and a transverse relative 
momentum squared p 2 . For the semileptonic rates 
the above mentioned relation takes the form 



dTsL(H Q ) 
dt 



dT SL 



(1) 



where dTg L /dt is the free quark differential decay 
rate, t = q 2 /ttIq, q being the 4-momentum of the W 
boson, and Rh q (t) incorporates the nonperturbative 
effects related to the Fermi motion of the heavy quark 
inside the hadron. The expression for dr® L /dt for 
the case of non-vanishing lepton masses is given e.g. 
in 0. R Hcj (t) in (|l|) is obtained by integrating the 
bound-state factor w(t, s) over the allowed region of 
the invariant hadronic mass Mx Q , (Q' = c or s) for 
the underlying transitions Q — » Q' : 



7 II Q ,{' 

R H Q {t) = J dsuj HQ (t,s), 



(2) 



where s 



uh q (t, s) 



M xJ m l and 



Trrrib 



Q n 

min[l,X2] 



dx\^(x,pf)\ 2 



(3) 



In Eq. (||) xo = wiq/Mhq , £1,2 = x q + /q ± , where 
Qo + |l| is defined in the Ht rest frame and 
qo ± |q| are defined in the Q quark rest frame. 
In Eq. (||) the region of integration is defined through 
the condition x\ < min[l, X2], i.e. s m in = p = 
(mQ,/mnf, s max = x Q 2 (1 - xq^/1) 2 . For other de- 
tails see 0. 

The structure of Eq. (jl]) suggests that in the 
limit of heavy quarks with infinite mass (i.e., tuq — > 
00 and ttiqi — > 00) one has 



dq uj(q 2 ,q ) = 1. 



(4) 



which means that the total inclusive width of the 
hadron is the same as the total inclusive width at the 
free quark level. The corrections to the free-quark de- 
cay picture are mainly due to the difference between 
the quark mass ttlq and the hadron mass Mh q as well 
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as to the primordial motion of the Q-quark inside 
Hq. These non-perturbative corrections vanish in 
the heavy-quark limit ttlq — > oo, but at finite values 
of the 6-quark mass a new parton description of in- 
clusive semileptonic decays, based on the constituent 
quark model, has been derived. Note that the quark 
mass rriQ appearing in|| is the constituent quark mass 
that differs from the pole quark mass mQ tPO i e that 
enters the OPE result for the inclusive widths. As a 
consequence, Eq.[j] generally contains the correction 
to first order in I/tuq expansion to the free quark 
decay width. In principle, it is always possible to re- 
define the constituent quark mass mq in such a way 
that the 1/rriQ correction is eliminated from the to- 
tal semileptonic width in accordance with the general 
statement of the OPE method. This is completely 
analogous to the way used to eliminate 1/mg cor- 
rections from the ACCMM model || and has been 
discussed within the considered approach in || for a 
particular case of the B — > X s j decay width. 

A priory, there is no connection between the 
equal-time (ET) wave function iu(kj) of a constituent 
quark model and LF wave function ip(xi,p? t ). The 
former depends on the center-of-mass momenta k^, 
while the latter depends on the LF variables Xi and 
Pit. However, there is a simple operational connec- 
tion between ET and LF wave functions |l0|. The 
idea is to find a mapping between the variables of 
the wave functions that will turn a normalized solu- 
tion of the ET equation of motion into a normalized 
solution of the different looking LF equation of mo- 
tion. That will allows us to convert the ET wave 
function, and all the labor behind it, into a usable 
LF wave function. This procedure amounts to a se- 
ries of reasonable (but naive) guesses about what the 
solution of a relativistic theory involving confining 
interactions might look like. 

We convert from ET to LF momenta by leav- 
ing the transverse momenta unchanged, pit = kit, 

and letting Xi = J=^-, where P + = Ylpt- ^ n 



the center-of-mass frame P + 
the free mass operator, Mq 



■Pi. 



and 



= ~(Xi 



P 2 t 



x t M§ 



-)M 



Mo, where Mq is 
^2 with cjj = 



(5) 



The momentum fractions and the transverse mo- 
menta obey the conservation law 



Now we obtain the LF wave function from 

d(k lz ,k 2zi k 3z ) 
tp{Xi,Pit) = -5-; r- • w(ki,k 2 ,k 3 ), 

d(xi,x 2 ,x 3 ) 

The straightforward calculation yields Jll| 

d(ki z ,k 2z ,hz) _ 
d(x 1 ,x 2 ,x 3 ) xix 2 x 3 M 

with 

W£ _ 1 

Xi ~ 2 



pi + m i 



xjMt 





It can be easily verified that 
1 

J [dxi] J [dpit\\4>{xi,p lt )\ 2 


- 1)*£>) = 1 

i i 

where 

[dxi] — dxidx2dxsS(^^Xi — 1), 

i 

[dpu] = d 2 p lt d 2 p2ttd 2 p 3t tS{^ pa) 

i 

provided 

ILd^Hkx^ka)! 2 
•<5(ki +k 2 + k 3 ) = 1. 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



Xi + X 2 + X 3 = 1, pit + p 2t + p 3 f = 



(6) 



Wave functions made kinematically relativistic in this 
fashion are used to calculate the form factors of 
heavy-to-heavy and heavy-to-light exclusive transi- 
tions Q. 

3. The Kalman d'Souza model 

As was mentioned above, we consider the 
problem of a constituent bound state formed by a 
heavy quark interacting with the lighter ones. Specif- 
ically, we consider the Kalman-D'Souza model for 
heavy baryons |l3] ]. 

The constituent quark masses are 

m u = md == 0.23 GeV, m s = 0.605 GeV 

m c = 1.961 GeV, m b = 5.637 GeV (13) 

We denote the light quarks (u,d, or s) by the labels 1,2 
and the heavy quark (c or b) by the label 3. For the 
cases of A and E baryons where two of three quarks 
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have equal masses we introduce the Jacobi relative 
coordinates 



/>= ^=(ri-r 2 ), A 



min + m 2 r 2 
mi + m 2 



r 3 (14) 



For the A c and baryons in the ground state, the 
predominant configurations (« 95%) of the heavy 
baryon wave functions contain no excitations in p and 
A. Hence 

3/2 



1 





n 






















exp 





„2 \ 2 



• exp 

where the factors a 7 (7 = p, A) are given by 



(15) 



then the momentum space wave function ip(h p ,]i\) 
normalized according to (0) is given by 



V>(k p ,k A ) = 



with 



1 



• exp 



3/2 



k 2 



1 



3/2 



k 2 



5* -P(-5*| (20) 



3«a 



(21) 



4. The heavy baryon light front wave func- 
tions 

We denote x=x 3 , y=xi, p3= pq, then 



MA 



with 
p p = 2 
and 

= 0.29 CeV. n 



mirri2 
mi + m 2 



3m 3 (mi 



2 (mi + m 2 + m 3 ) 



(0) 
A 



0.32 GeV. 



(16) 



(17) 



(18) 



In Table 1, we compare the theoretical predictions 
of the Kalman-D'Souza model with the experimental 
data. 

Table 1. The theoretical and experimental 
masses of the heavy baryons and the parameters Xq, 
a p , and a x used in the calculations. The theoretical 
masses are taken from fllS], while the experimental 



ones are taken from ]14fl (for H c and Ab) and from 
PI for ~ b 



Hq 


M th 




x 


a p 


a\ 


A c 


2.280 


2.285 ±0.0006 


0.86 


0.29 


0.40 


S c 


2.454 


2.453 ±0.0006 


0.80 


0.29 


0.40 




2.461 


2.466 ±0.0014 


0.80 


0.29 


0.40 


A b 


5.640 


5.624 ±0.009 


1.00 


0.29 


0.41 


s 6 


5.960 




0.95 


0.29 


0.41 


2b 


6.045 




0.93 


0.29 


0.41 



l<p(z, Pt)l : 



dy / dp lt 



• yM 

• (m 4 

where 



kl = ( Pi 



Pq* + m Q 



Mn 



(Pot 



(1 - x - y y M 



pit + m i 

y 2 M 



l^(k p ,k A ) 



mi + 7712 

1 



PQ 



K = Pq = Pt + t xM 



Pt +m 2 Q 
xM 



(22) 



(23) 



(24) 



In Fig.l we show the two dimensional plot 
of the distribution function n (x,p 2 ). In Fig. 2 are 
shown the one dimensional distribution functions 
F Hq (x), x=x 3 . 



Fh q (x) = 



1 3 
dxi J dx 2 J JJdpi, 



•\i>H Q (xi,pu)\ 2 

.sfcPitWY,^- 1 ) 



(25) 



We now introduce the charge conjugate mo- 



menta 



k P = i(ki - k 2 ), k A = "^i+k2)-2m 9 k 3 i (ig) 
2 Zrriq + tuq 



for H=A,S and Q=b,c. 

As illustrated in Fig. 2 the distributions fc c and fe^are 
modified by the presence of the s quark instead of 
u(d) one. In Fig. 3 we show the non-perturbative 
factor Rfj Q (t) for Ab,Sb, A c ,and 5 c baryons. 



Figure 3. The bound-state factor Rh q in Eq. (H). 



Figure 1. Two-dimensional distribution function 
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Figure 2. One-dimensional distribution functions 
Fjj (x) for A c and A baryons. 



5. The results 

We have evaluated Eqs.(Q-||) in case of the 
inclusive semileptonic decays Hq — > Hqi£i>i adopt- 
ing the heavy quark distribution function Fh q (x) 
from ( p5| ) with the Gaussian ansats of the three- 
quark wave function (|l5|), that corresponds to the 
constituent quark model of ref . |l3[ . The parameters 
in (15) have been taken from Table 1. Our results 



for the SL decay widths of H c and H b baryons and 
electron branching ratios TsL(HQ)/T tot (HQ), where 
T tot (H Q ) = 1/t(H q ) is the total width of H Q baryon, 
are collected in Table 2. For the lifetimes we take the 
values from the recent PDG publication Il4j r(A c ) = 
0.206 ± 0.012 ps-\ r(H+) = 0.33 Ig^ ps~\ 
r(A b ) = 1.229 ± 0.08 ps" 1 , r(S 6 ) = 1.39 ± 0.3 ps _1 . 

For the c-flavored baryons, there exists only 
one inclusive value for the semileptonic branching ra- 
tio BR(A C — > X s ev) = 4.5 ± 1.7%, the corresponding 
exclusive semileptonic branching BR(A C — * Kev) be- 
ing 2.0 ±0.6% fL§. For S+ -> E°ev there exists only 
exclusive branching ratio 2.3~^'l%. For the 6-flavored 
baryons, there exists the semileptonic branching ra- 
tio BR(A b -> A c £v e + anything) = 7.9 ± 1.9% which 
is dominated by the charm baryons. In Table 2, in 
parentheses are also shown the results obtained us- 
ing OPE with the pole quark masses m c = 1.6 GeV, 
m b = 4.8 GeV and \V cb \ = 0.039 [@. We should 
mention a rather good agreement with the OPE pre- 
dictions for the semileptonic rates in spite of the large 
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difference of the constituent and pole quark masses 
that produces the large difference in dTg L /dt ~ rriQ. 
This difference is compensated by a factor R(q 2 ) in 
Eq. (|). 



Table 2. The predicted value of the inclusive 
semileptonic widths (in units of ps _1 ) and the cor- 
responding branching ratios in %. The experimental 
data are from pi. \V cb \ = 0.04, \V SC \ = 1.00. 



Hq 


T{X Q ,ev e ) 


T{X q ,tv t ) 


BR{X Q ,ev e ) 


Ac 


0.267 (0.233) 




5.50 


1 — 'C 


0.236 (0.233) 




7.79 


A b 


0.086 (0.070) 


0.019 


10.57 




0.092 (0.070) 


0.021 


12.79 
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6. Conclusions 

This is the most consistent quark model to- 
date describing baryons in the sense that a single set 
of parameters is used for the whole spectra and then 
subjecting all of their properties namely the masses, 
radii and their wavefunction structures to a strict test 
by a decay model. The spectroscopy is successfully 
calculated using far less parameters than any other 



model 13 1. Decays of baryons are calculated [|5j. Now 
in this paper we have successfully computed the in- 
clusive semileptonic decays of the heavy baryons Hq 
by incorporating the LF formalism of 
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